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Structure and reaction studies with a method of antisymmetrized molecular dynamics 
(AMD) were reviewed. Applications of time-independent and time-dependent versions of 
the AMD were described. In applications of time-independent AMD to nuclear structure 
studies, structures of neutron-rich nuclei such as Be, C, Ne, and Mg isotopes were described 
focusing on cluster aspects. Important roles of valence neutrons were discussed. The results 
suggested a variety of cluster structures appear also in unstable nuclei as well as in stable 
nuclei. Deformation and cluster phenomena \n Z N nuclei in p- and sd- shell regions 
were also discussed. Applications of time-dependent AMD contain various topics such as 
fragmentation in heavy-ion collisions as well as nuclear responses. The AMD calculations 
successfully describe multifragmentation which is one of the remarkable phenomena in heavy- 
ion collisions. The approach is able to link reactions and nuclear matter properties. The 
studies suggested the important balance between single-particle motions and correlations to 
form clusters and fragments. 

Contents 



1. Introduction [2] 

2. Antisymmetrized molecular dynamics (AMD) [3] 

2.1. Basic formulation of AMD d 

2.2. Basis AMD and its Extensions in applications [71 

2.3. Branching in time evolution [9] 

3. Applications of time-independent AMD method to nuclear struc- 
ture [11] 

3.1. Molecular structures in Be and Ne isotopes [TT] 

3.2. Three-body cluster states in ^^c, i^B, and ^^C [H 

3.3. Spectroscopy and exotic phenomena in the island of inversion [18] 

3.4. Superdeformation in sd-shell nuclei [23] 

3.5. Decoupling between Proton and neutron deformations [27] 

4. Applications of time-dependent AMD method to nuclear response 
and reaction [29] 

4.1. Dipole resonances [30] 

4.2. Radial oscillations [31] 

4.3. Multifragmentation in expanding systems [33] 

4.4. Fragmentation in collisions of light nuclei [36] 

4.5. Statistical properties of excited systems [33 

4.6. Symmetry energy effects in heavy-ion collisions [40] 



typeset using -pTT^TeX.cIs (Ver.0.9) 



2 



Y. Kanada-En^yo, M. Kimura, and A. Ono 



5. Summary and perspective 1421 

§1. Introduction 

A nucleus is a finite quantum many-body system consisting of protons and neu- 
trons interacting via nuclear forces. Its ground state has shell structure, in which 
nucleons move almost independently in an averaged field (mean field) analogously 
to an atomic system. The shell structure and excitation modes associated with 
single-particle and collective motions are important facets of nuclear system. 

On the other hand, in nuclear systems, one may also find unique characters 
different from atomic systems. One of the important differences is that a nucleus is 
a self-bound system formed by attractive nuclear forces. Because of the attraction, 
spatial correlations among nucleons can be rather strong, and therefore, assembling 
and disassembling of nucleons occur in various ways. Furthermore, the saturation 
property of nuclear system, where binding energy per nucleon and central density 
are almost constant independently to the mass number, implies that assembling and 
disassembling of nucleons can take place with a small excitation energy. As well 
known, this characteristic of nuclear system manifests itself as cluster structures in 
which a nucleus is divided into several subunits (clusters) and nucleons are confined 
within each cluster. The important roles of assembling and disassembling of nucleons 
continue to higher excitation energies as in intermediate-energy heavy-ion collisions 
where a lot of clusters and fragment nuclei are produced from a hot source whose 
excitation energy is typically comparable to the binding energy of a nucleus. 

In spite of importance of cluster aspect in nuclear systems, usual mean-field 
approaches often fail to describe those cluster phenomena because they take into ac- 
count insufficiently many-body correlations, which are essential in cluster formation. 
The theoretical method of antisymmetrized molecular dynamics (AMD)[ll>'§'ll}'EI'[l}'§'[Zl 
has been proposed in studies of heavy-ion collisions. AMD describes nuclear many- 
body systems by antisymmetrized products of Gaussian wave packets of nucleons 
and incorporates quantum effects and nucleon-nucleon collisions. It has described 
fragment formation in heavy-ion collisions successfully. AMD has been proved to be 
a powerful tool also for nuclear structure study. All centers of the Gaussian packets 
are independently treated as variational parameters in the AMD framework, and it is 
possible to describe various cluster structure wave functions as well as independent- 
particle motion in a mean field without a priori assumption. Thus, AMD has been 
applied to investigate various phenomena in nuclear structure and reactions. 

In the early days, AMD studies were hmited to light systems. This is because 
computational cost increases rapidly in proportion to A^^^ {A is the mass number) 
due to the non-orthogonality of single-particle wave functions. However, AMD cal- 
culations have developed remarkably toward various nuclear systems owing to the 
rapid progress of computational facilities. For instance, it enabled the reaction stud- 
ies by AMD up to Au+Au collisions {A ^ 400) |^ the structure studies up to p/-shell 
nuclei {A ^ 40)Hi and variational calculation after the angular momentum projection 
that covers up to very high excitation energy© Further high-performance comput- 
ing will extend the subjects of AMD studies and enable more sophisticated AMD 



Antisymmetrized molecular dynamics 



3 



calculations. 

AMD studies have revealed that cluster phenomena emerge widely in various 
nuclear systems. Many exotic and novel features of clustering have been discovered 
and the concept of the clustering has been renewed and extending now. Today, it 
is well established that clustering is an essential aspect of nuclear many-body sys- 
tems as well as the mean- field aspect. Coexistence of cluster and mean- field aspects 
brings out rich phenomena to nuclear many-body systems as functions of excitation 
energy and isospin degrees of freedom (Fig. [T]). As the excitation energy increases, 
one may see transitions from mean-field to cluster structures. In deeply bound sys- 
tems such as low-lying states of stable nuclei, mean-field effects are rather strong. 
However, even if a nucleus has a shell-model-like structure in its ground state, devel- 
oped cluster structures appear in excited states near the corresponding cluster-decay 
threshold energy (so-called Ikeda's threshold rule^^^). Above the threshold energy, 
further remarkable cluster phenomena such as alpha decays, molecular resonances, 
and fission etc. are known. The coexistence and competition between the clustering 
and mean field have been studied by AMD in a unified way within a single theoret- 
ical framework. When the excitation energy increases further, nuclear systems may 
enter the region of nuclear liquid-gas phase transition. AMD studies have confirmed 
the phase transition by obtaining caloric curves for equilibrium systems. The link 
between the phase transition and copious fragment formation in heavy-ion collisions 
is now clearer with the unified description by time-dependent AMD. Another im- 
portant degree of freedom in nuclear systems is the isospin asymmetry, i.e., neutron 
(proton) excess, which is one of the major directions in recent nuclear physics. In 
the neutron-rich and proton-rich domains, the saturation law of energy and den- 
sity is broken. Because of the unbalanced proton-neutron ratio, we may encounter 
many novel cluster phenomena in isospin asymmetric systems. Indeed, such exotic 
clustering phenomena as cluster structures in neutron-rich Be isotopes and isospin 
fractionation/distillation at liquid-gas separation in fragmentation reactions have 
been investigated with AMD. 

In this paper, we review the AMD approach and its applications to nuclear 
structures and reactions. In the next section, the formulation of AMD is described. 
Applications of the time-independent version of AMD to static problems of nuclear 
structures are explained in ^ and those of the time-dependent version to dynamical 
phenomena such as nuclear responses and nuclear reactions are described in Q 
Finally, a summary and perspectives are given in ^ 

§2. Antisymmetrized molecular dynamics (AMD) 

As mentioned above, cluster aspect is one of the essential features of nuclear sys- 
tems as well as mean- field aspect. Originating in coexistence of two kinds of nature, 
cluster and mean-field aspects, a variety of phenomena arise in nuclear many-body 
systems. To investigate rich phenomena concerning cluster aspects such as cluster 
structures and multifragmentations, a theoretical framework that can describe both 
the cluster and mean-field features is required. The AMD method has been proved 
to be one of the powerful approaches in description of those features. In this section. 
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Fig. 1. Schematic figure for rich phenomena in nuclear systems. 



we briefly review formulation of AMD. After describing basic formalism of AMD, we 
explain extended methods of AMD which have been applied to nuclear structure and 
nuclear reaction studies. For detailed formulation of the AMD method, the reader 
is refereed to Refs. Q,[5|,[6]),[7|. 

2.1. Basic formulation of AMD 
2.1.1. AMD wave function 

In the AMD framework, an basis wave function for an A-nucleon system is 
expressed by a Slater determinant of Gaussian wave packets; 

1 



^amd(^) 



(2-1) 



where the ith. single-particle wave function is written by a product of spatial (0), 
intrinsic spin (%), and isospin (r) wave functions as follows. 



(2-2) 
(2-3) 
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(2-4) 



The spatial part (j)Zi of the ith single-particle wave function is represented by a 
complex variational parameter Zicr with a = x^y^z which indicates the center of 
the Gaussian wave packet. The spin part Xi is parametrized by a complex number 
parameter, The isospin function is fixed to be up (proton) or down (neutron). 
The width parameter z/ takes a common value for all nucleons. It is chosen to be 
an optimum value for the studied system. Accordingly, an AMD wave function is 
expressed by a set of variational parameters, Z = {Zi, Z2, • • • , Z^,^i,^25 • • • ^Ca}- 
These parameters indicate centers of localized Gaussians and spin orientations, which 
are treated independently for all nucleons. That is to say, a system written by a single 
AMD wave function is specified by the configuration of single-nucleon wave packets 
in the phase space and their spin orientations. In the simplest version of AMD, the 
spin part Xi is sometimes fixed to be up or down, and only parameters Zi, Z2, • • • , Z^ 
for Gaussian centers are treated as variational parameters. In time-dependent AMD, 
intrinsic spin orientations are usually fixed. In the case of fixed spins, we redefine 
Z = {Zi, Z2, • • • , Z^} by omitting the spin labels 

In the AMD wave function, all single nucleons are treated independently as 
localized Gaussians. Although any constituent clusters are not assumed a priori^ 
multi-cluster structures can be described by grouping of single-nucleon Gaussian 
wave packets in the spatial configuration. On the other hand, if all the Gaussian 
centers gather around a certain position, the AMD wave function becomes equivalent 
to a harmonic oscillator shell-model wave function around the position due to the 
effect of antisymmetrization. Thus, the model space of AMD can describe both 
cluster and mean-field features with assembling and disassembling of Gaussian wave 
packets. If a system favors a specific cluster channel, such a cluster structure will be 
automatically obtained in energy variation or in dynamics. 

The AMD wave function is quite similar to a wave function of fermionic molec- 
ular dynamics (FMD)"^'[llf where more generalized wave functions are adopted. In 
applications of AMD to nuclear structure and reaction studies, the description has 
been improved by superposing many AMD wave functions or by introducing stochas- 
tic processes rather than by choosing more general single-particle wave functions. 

2.1.2. Equation of motion 

In the time- dependent versions of the AMD method which have been applied to 
dynamics of nuclear systems, the time evolution of the variational parameters Z are 
determined by the time-dependent variational principle. The equation of motion for 
Z derived from the time-dependent variational principle is 




(2-5) 



where a, p = x^y^z are the labels for the components of Z^ (z = 1, 2, . . . A). When 
the spin wave function is a variable, it may be regarded as the fourth component 
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Zi4 = ^i- The expectation value of the Hamiltonian H is given by 

n,rr, _ (^AMD(^)lffl^AMD(^)) 

(^amd(^)I^amd(^)) • ^ ' 

A positive definite Hermitian matrix 

Ciajp ^ gJlz.^ H<Pamb{Z)\^amd{Z)) (2-7) 

appears in the equation of motion, suggesting that the variables Z are not canonical 
coordinates. 

2.1.3. Energy variation 

To get an optimum solution for the energy minimum state, the energy variation 
is performed. Namely, the variational parameters Z are optimized to minimize the 
expectation value of the Hamiltonian in the AMD model space. We introduce the 
following frictional cooling equation. 



The parameter A is an arbitrary real number and ii is an arbitrary negative real 
number. It is easily proved that the energy of the system decreases as time develops 
due to the frictional term iji. In the energy variation, the matrix Ci^jp can be 
replaced with SijS^p and then, the frictional cooling method with A = becomes 
equivalent to the steepest decent method. In both cases of Ci^jp^ we obtain the 
optimum set of parameters Z that gives the AMD wave function for the minimum 
energy state in the model space after enough cooling time (iteration steps). 

2.1.4. Hamiltonian 

The Hamiltonian for an A-nucleon system consists of the kinetic energy, the 
nuclear and Coulomb force terms, 

^ = E + E + E - (2-9) 

Here ti is the kinetic energy. The energy of the center-of-mass motion Tg is subtracted 
exactly because the total wave function can be separated into the internal wave 
function and the center-of-mass wave function. In applications to heavy ion collisions, 
spurious kinetic energy of the zero-point oscillation of fragment mass centers is also 
subtracted from the Hamiltonian. For the effective two-body nuclear force Vij^ finite- 
range forces such as forces supplemented by GSRS form spin-orbit forces}i2f 
Gogn}J^'^ forces, and forces are used. Finite-range two-body forces 
with zero-range three-body forces such as the modified Volkov forces^^' are also used. 
The Coulomb force v99^^^^^ is approximated by a sum of seven Gaussians. 

These effective forces are phenomenological ones constructed to describe low- 
energy properties of nuclear structure. In heavy-ion reaction, residual interactions 
contribute to nucleon-nucleon collisions which are incorporated by stochastic collision 
process in the AMD framework as explained later. 
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2.2. Basis AMD and its Extensions in applications 

In the early stage of AMD studies, simple versions of the AMD method have 
been applied to reaction and structure studies j3 'HI* '0 and later, the AMD method has 
been developed to many extended versions (Refs. [4]),[6]),[7|) and references therein). 
In this section, we explain the formulation of the basic AMD method and some 
advanced versions for structure study and those for reaction study. 

An basis AMD wave function is given by a single Slater determinant. Gen- 
erally, many-body wave functions for quantum systems should be expressed by a 
superposition of many Slater determinants. Restriction of the model space within 
a single Slater determinant is the limit of a mean-field approximation. To incor- 
porate beyond-mean-field effects, superposition of Slater determinants is essential. 
Firstly, in structure study, parity and angular-momentum projections, which are 
done by superposition of Slater determinants, are essential to describe properties of 
energy-eigen states. In second, superposition of Slater determinants is significant to 
improve wave functions by taking into account quantum fluctuation around a mean- 
field, many-body correlations, and spin-parity projections. It is also necessary to 
describe excited states orthogonal to lower states. 

In applications of the AMD method to structure study, the projections and 
superposition of AMD wave functions are practically performed. In applications to 
heavy-ion reactions, however, emergence of multiple reaction channels is introduced 
by stochastic branching processes, neglecting quantum interference between different 
channels each of which is described by an AMD wave function. 

2.2.1. Projections and superposition of AMD wave functions 
The parity-projected AMD wave function is given as 

I^amd) ^ ^'^I^amd(^)) = ^^I^amd(^)), (2-10) 

where = 1 ± is the parity projection operator. The angular- momentum 
projected AMD wave function is written as 

\Kk) = Pmk\^amd{Z)) = J df2Di^K{f2)R{f2)\0AMD{Z)). (2-11) 

Here Dj^^{Q) is the Wigner's D function and R{Q) is a rotation operator with 
respect to Euler angle i7. As clearly shown, the angular-momentum projected state 
is expressed by a linear combination of wave functions rotated from the intrinsic 
AMD wave function ^amd(^) with the weight function, D'^j^. The matrix element 
of a tensor operator T^, where k is the rank and q is the z-component, can be 
calculated to be 

{Pi,K^AMD{Z)\fg'\Pii,K>^AMD{Z')) = ^^{J' M' kq\J M) 

xJ2{J'f^kiy\JK) I dQD-i;^,{Q){$AMD{Z)\f^R{Q)\$AM^{Z')). (2-12) 
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In practical calculations, the integrations with respect to f2 is performed by numerical 
integration on grid points of angles f2 = (^1,^25^3)- In calculations of expectation 
values for observable operators such as Hamiltonian, radii, moments, and transitions, 
AMD wave functions are projected to parity and angular-momentum eigenstates. 
In the usual AMD calculations for structure study, the parity projection is done 
before energy variation while the angular- momentum projection is performed after 
the energy variation, i.e., variation before projection (VBP). 

Superposition of independent AMD wave functions is useful to improve wave 
functions and it is essential in description of excited states to satisfy orthogonality 
between energy levels. Let us consider superposition of independent AMD wave 
functions ^amd{Z^^^) (A: = 1, • • • , fcmax) (^max is the number of adopted basis AMD 
wave functions). Superposed wave functions for states is written as 

K^) = E^-^^^^I^mx^amd(^('^)), (2-13) 

kK 

where = Pmk^^- H^^'^ the values for the coefficients c^^kJK are determined 

by the variational principle, 

<5{«±|^|<±) - 6„«±|<±)} = 0, (2-14) 

which is equivalent to diagonalization of the norm and the Hamiltonian matrices 



and leads to the Hill- Wheeler equations. The K sum in Eq. (2-13) stands for the 
i^^- mixing. 

To adopt efficient AMD wave functions as basis wave functions for the superpo- 
sition, constraint methods, which are often used in a generator coordinates method 
(GCM)jl3'[ll'lll}@ or stochastic variational methods^^^ are apphed in the AMD 
framework. Here we consider a constraint {g) — g. By changing the constraint value 
as ^ = gi,g2r'' ^9kma^i the minimum energy state ^amd(^) the AMD model 
space is obtained for each constraint value by the constraint energy variation after 
parity projection. Then the obtained wave functions ^amd(^^) ~ • • • 5 ^max) are 
superposed, and coefficients are determined by the diagonalization. This method cor- 
responds to a GCM calculation for a generator coordinate g when an enough number 
of the basis wave functions for different values gk are taken into account. This method 
(called AMD+GCM) is useful, in particular, for study of excited states. For the con- 
straints, the nuclear quadrupole deformation parameters are often used Hll'llll* '113 In 
the AMD+GCM method, variation is done before the angular-momentum projec- 
tion. 

It is also efficient to perform variation after the angular-momentum projection 
as well as the parity projection, especially, for study of excited states. Namely, the 
energy expectation value for a parity and angular- momentum projected AMD wave 
function, 

^ ^ {P^1,<pAMB{Z)\H\P'^i,<PAMB{Z)) 
{PMK^AMD{Z)\PitK^AMD{Z)) ' 

is minimized with respect to the variational parameters Z by using the frictional 
cooling method. Firstly, for a given spin and a parity J^, the wave function for 
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the lowest state is obtained by the energy variation. Then, a wave function 
for a higher state (J^) is provided by varying Z to minimize the energy for the 
component of an AMD wave function orthogonal to the lower states (J^, • • • , J^_i) 
which are already obtained. In the present paper, we call the variation after spin- 
parity projection "VAP".f^ 

2.2.2. Extensions for single-particle wave functions 

In a basis AMD wave function, a single-particle wave function is written by a 
spherical Gaussian. In extended versions of AMD, a single-particle wave function 
is written by a deformed Gaussian or superposition of different-range Gaussians 
to improve single-particle wave functions.^^''^^''^^^''^ In the method of deformed- 
basis AMD,^^^'^^' triaxially deformed Gaussians are employed as single-particle wave 
packets instead of spherical ones. 



f>z,(rj) ocexpj- ^ M^j^ - 



a=x,y,z 



(2-16) 



where the width parameters, i^y and z/^, take different values for each direction. 
They are determined in the energy variation to optimize energy of a system. By using 
this deformed basis, it is possible to successfully describe coexistence (or mixing) of 
cluster and deformed mean- field structures, which are essential especially in heavy 
systems. 

Superposing different range Gaussians is another method to improve single- 
particle wave functions of an AMD wave function as is done in the FMD method H^f'lHZf 

Instead of adopting deformed or superposing Gaussians, a method of stochas- 
tic branching on wave packets is used to describe the diffusion and deformation of 
single-particle motions in the time-dependent version of AMD for nuclear reaction 
calculationiU'lllI® as explained in the next subsection. 

2.3. Branching in time evolution 



Initial State 




Fig. 2. A schematic picture of the quantum branching processes for multichannel reactions. 



The description of dynamics of excited nuclear many-body systems, such as in 
heavy-ion collisions, is a highly quantum-mechanical many-body problem. If the 
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many-body time-dependent Schrodinger equation is solved for an initial state which 
may be roughly approximated by a single Slater determinant, the intermediate and 
final states will be a very complicated states containing a huge number of reaction 
channels corresponding to different fragmentation configurations, as illustrated in 
Fig. [2j The AMD model for reactions respects the existence of channels, while it 
neglects some of the interference among them. Namely, the total many-body wave 
function \^{t)) is approximated by a mixed state 

where each component is represented by an AMD wave function \^{Z)) with a time- 
dependent weight w{Z^t). 

Interference between quite different components is not of practical importance 
because the matrix elements are negligible for the Hamiltonian and other operators 
of our usual interest. Furthermore, when we adopt mean field approximation in 
some way, we need to take care of the spurious nonlinearity introduced by the ap- 
proximation. A many-body state composed of different reaction channels cannot be 
described by a single mean field. Our strategy is to suitably decomposed the state 
in such a way that the mean field approximation is valid in each component. 

Let us consider the motion of a nucleon in the system, ignoring the Pauli prin- 
ciple for the moment. In mean field approximation, the one-body density matrix is 
given by a pure single-particle wave function p{t) = \ilj{t)){i(j{t)\ at any time. Un- 
der the mean field, the phase space distribution will spread in some directions and 
may shrink in other directions. In contrast, in the real time evolution, the condi- 
tion = p{t) should no longer hold due to many-body correlations. Namely, 
decoherence should occur on the single-particle state so that it turns into a mixed 
state. In the AMD approach, decoherence is taken into account by splitting the 
wave packet in such a way that the spreading of the distribution in the mean field is 
respected while the shrinking is discarded. By the decoherence of the single-particle 
states, the decomposition (or branching) of the many-body state is induced natu- 
rally. This approach of using compact wave packets is advantageous for the case with 
many fragmentation channels because it is free from spurious coupling of different 
fragmentation channels. 



Instead of directly considering the weight function w{Z,t) in Eq. (2-17), we 
solve a stochastic equation of motion for the wave packet centroids which may be 
symbolically written as 

j^Zi = {Zi, n}pB + (NN coll) + AZi{t) + II (Zi, n'). (2-18) 
The first term {Z^,?-^}pb is that in the deterministic equation of motion [Eq. 



(2-5)] derived from the time-dependent variational principle. 

The second term represents the effect of stochastic two-nucleon collisions, where 
a parametrization of the energy-dependent in- medium cross section is adopted. The 
collisions are performed with the "physical nucleon coordinates" that take account 
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of the antisymmetrization effects, and then the Pauh blocking in the final state is 
automatically introducedH'H 

The third term AZi{t) is a stochastic term for the wave packet splitting men- 
tioned abovellf'lllf'llHl The change of the width and shape of each wave packet is 
calculated by solving the Vlasov equation (for some time period) with the same ef- 
fective interaction as for the term {Z^,'H}pb. An essential ingredient here is how 
long the coherent single-particle motion is solved before decoherence. The properties 
of fluctuations AZi{t) are determined depending on this time scale called coherence 
time T. When the wave packet splitting was first introduced into AMD^^'^i^ the 
limit of T ^ was taken, for which the decoherence effect is maximum. The choice 



of a finite coherence time has been formulated in Ref. 29), and a reasonable choice 
may be to consider decoherence for a nucleon when it collides with another nucleon. 
This choice is symbolically denoted by r = tnn, though r is not a constant number. 

The last term ji{7ii^l-L') is a dissipation term related to the fluctuation term 
AZi{t). This term is necessary in order to restore the energy conservation that is 
violated by the fluctuation term. The coefficient fi is given by the condition of energy 
conservation. However, the form of this term is somehow arbitrary. The variables Z 
are shifted to the direction of the gradient of the energy expectation value T-L under 
the constraints of nine conserved quantities (the center-of-mass variables and the 
total angular momentum). When the number of other nucleons around the nucleon 
i within a phase space radius of |Zj — Z^| < 2.5 is more than a certain number 
NcQ, other twelve global one-body quantities (monopole and quadrupole moments in 
coordinate and momentum spaces) are also included as constrained quantities for the 
dissipation corresponding to the fluctuation AZi{t). The latter constraints should 
be imposed because the one-body time evolution has already been considered by 
{Zi,T-L}pB + AZi{t). Empirically, A^c9 has been chosen between 5 and 15. 

A summary of the complete formulation of AMD for reactions can be found in 
Ref.|4]). 

§3. Applications of time-independent AMD method to nuclear 

structure 



We here discuss some topics investigated with the time-independent AMD method 
focusing on cluster aspects. 

3.1. Molecular structures in Be and Ne isotopes 

Cluster structure of Be isotopes is one of the fascinating subjects of unstable 
nuclei. A 2a-cluster core is favored in neutron-rich Be isotopes as well as ^Be whose 
ground state is a 2a resonance state. The low-lying states of neutron-rich Be isotopes 
are described well by a molecular-orbital picture based on a 2a core and valence 
neutrons moving around the 2all2l'[13@@'[lll'l^ In contrast to the molecular- 
orbital structures in low-lying states, developed di-cluster states such as ^He+^He 
in -'^^Be have been suggested in highly excited states7^'^^''^^''^^^'^^^ There, valence 
neutrons are moving not around the whole system but around one of two a clusters. 
It means that a variety of cluster structures coexist in neutron-rich Be isotopes where 
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valence neutrons play important roles. 

The molecular-orbital picture has been extended also to Ne isotopes such as 
^^Ne and ^^Ne based on an ^^O+a-cluster core and valence neutrons in molecu- 
lar orbitalsI3@'E2l Di-cluster states like ^^0+a-cluster states in ^^Ne is another 
attractive subject Fl 

3.1.1. Molecular-orbital structure 

The idea of the molecular orbitals surrounding a 2a core was suggested in ^Be 

with a 2a+n cluster model^^^ in 1970's. In 1980's and 1990's, molecular-orbital 

models were applied to neutron-rich Be isotopes and succeeded to describe rotational 
bandsi2J'3i|@,[33| 

In a 2a system, molecular orbitals are formed by a linear combination of p 
orbits around two a clusters. In neutron-rich Be isotopes, valence neutrons occupy 
the molecular orbitals around the 2a core. The negative-parity orbital is called V 
orbital', while the longitudinal orbital with positive parity is V orbital' (Figjs]). bmce 
the a orbital has two nodes along the a-a direction, it gains the kinetic energy as 
the 2a cluster develops. The energy gain of the a orbital in the developed 2a system 
results in the intruder configurations of the ^^Be and ^^Be ground states. In other 
words, it is the origin of the breaking of the neutron magic number = 8 in the Be 
isotopes. 

In analogy to neutron-rich Be isotopes, molecular-orbital structures in Ne iso- 
topes have been suggested from the experimental systematics.*^ Indeed, AMD 
calculation has predicted the presence of the molecular-orbital bands with 
cluster core surrounded by two valence neutrons in the a orbital0® Different from 
Be isotopes, a orbital is a linear combination of sd orbits around ^^O and p orbits 
around a and it results in a p/-shell like molecular orbital. Another difference is the 
parity asymmetry of the core, that produces parity doublet of the molecular bands. 
Similar molecular-orbital structures have been also suggested for F isotopes© De- 
tailed discussions are given in later sections. 

3.1.2. Cluster structures in neutron-rich Be isotopes 

Cluster structures of Be isotopes have been intensively investigated in many theo- 
retical works with cluster models,^^''^^^'®© molecular-orbital models, ^^''^■'^''^^''^^''^ 
and AMdJ§®@®'I13 In the cluster and molecular-orbital models, the existence 
of two a clusters are a priori assumed. On the other hand, AMD does not rely 
on model assumptions of the existence of cluster cores. Nevertheless, the results of 
AMD calculations indeed indicate the appearance of the 2a core surrounded by va- 
lence neutrons in the molecular orbitals in low-lying states of Be isotopes. It means 
that the formation of the 2a core and molecular orbitals has been theoretically con- 
firmed by those AMD calculations without assuming clusters. Here we discuss the 
cluster structures of neutron-rich Be isotopes based on AMD calculations.^ 'Ell 

The systematic study of the ground and excited states of Be isotopes was per- 
formed with VAP calculations in the AMD model (AMD-VAP). Many rotational 
bands having the 2a core structure were obtained in the theoretical results. In 
Fig. [5) density distributions of the intrinsic wave functions for the band-head states 
of ^^Be and ^^Be are shown. As is seen, the proton-density distribution indicates 
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the formation of the 2a core, while the neutron density distribution exhibits the 
behavior of valence neutrons around the 2a core. 

In analysis of single-particle wave functions in the AMD wave functions, valence 
neutron orbits in low-lying states of Be isotopes were found to be associated with 
the molecular orbit als around the 2a core. In Fig. |4l we show schematic figures of 
the cluster states suggested in ^^Be, ^^Be, and ^^Be. In the figure, we show the 
number of neutrons occupying the cr-like orbit als, which have dominant positive- 
parity components and are regarded as the a orbital. The experimental value for 
the excitation energies of the corresponding states are also shown in the figure. 

In -'^^Be, the valence neutron configurations of ■'^^Be(O^), ■'^^Be(l~), and ■'^^Be(02") 
can be regarded as tt^, ttct, and configurations meaning two neutrons in the tt 
orbitals, one neutron in the tt and the other neutron in the a, and two neutrons 
in the a, respectively. Spatial distributions of the single-particle orbits of valence 
neutrons in -^^Be(O^) and ■'^^Be(02") shown in Fig. [sj^a), where the vr-like or- 
bital and the a-like orbital are clearly seen. Similarly to ^^Be, it was found that 
the ^^Be(l/2+),iiBe(l/2-), and ^^Be(3/22 ) states correspond to ttV, tt^ and Tva^ 
configurations, while the ■'^^Be(O^), ■'^^Be(02"), ^'^^^{^i) states are roughly inter- 
preted as TT^cr^, 7T^ and n^a^ configurations. 

Interestingly, the degree of the 2Q:-cluster development strongly correlates with 
the number of valence neutrons in the a orbital. Namely, the 2a cluster develops as 
the neutron number in the a orbital increases. It is easily understood because the 
single-particle energy of the a orbital decreases because of the kinetic energy gain 
in largely distant 2a systems. The enhancement of the 2a cluster with neutrons in 
the a orbital is consistent with the arguments in Refs. 31), 32), 33). On the other 
hand, as the neutron number in the tt orbitals increases, the cluster structure tends 
to weaken. 

Another interesting characteristic in Be isotopes is the breaking of neutron 
magicity in ^^Be and ^^Be. The breaking of the p shell for the neutron magic num- 
ber = 8 in ^^Be has been experimentally known from the unnatural parity 1/2^ 
ground state, while that in ^^Be has been suggested from slow /3 decay. Those 
exotic features of ^^Be and ^^Be can be understood from the molecular orbital pic- 
ture. The ground states of ^^Be and ^^Be are considered to have dominant intruder 
configurations with a-orbital neutron(s) instead of normal Ohuj configurations. The 
ground 1/2+ state of -'^-'^Be corresponds to the tt^ct configuration, while -'^^Be(O^) is 
the intruder state Tv^a^ in terms of molecular orbitals. In the one-center shell-model 
limit, the tt and the a orbitals correspond to the p and sd orbits. Therefore, in the 
ground states, -'^-'^Be(l/2^) and -'^^Be(O^) have dominant Ihuj and 2/icj configurations, 
respectively, indicating the vanishing of the N = 8 magic number in ^^Be and ^^Be. 
The breaking of the neutron shell in neutron-rich Be isotopes is caused by the lower- 



ing a orbital in the developed 2a structures as discussed in Refs. j6|),38),47). Again, 
the a orbital in the 2a structure plays an important role. 

In addition to molecular-orbital structures in such low-lying states, the AMD 
results for ^^Be suggested molecular resonant states having di-cluster ^He+^He and 
^He+o; structures in highly excited states. The result is consistent with the ex- 
perimental observations of cluster states in He+He break-up react ion^^@ and also 
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with theoretical suggestions by cluster model calculations HZl'IMI '113 




Fig. 3. (a) Density distributions of the single-particle wave functions for valence neutrons in 
^°Be(0+) and ^°Be(0+).^^' Schematic fi gure of the molecular orbit als, tt and a orbitals around 
the 2a core are also shown at the bottom, (b) Density distribution of the excited band {K^ = 0~) 
in ^^Ne obtained by AMD.^°' The middle and top figures show the density distribution of the 
single-neutron wave function of the highest single-particle level. The matter density of the total 
system is displayed at the bottom of the box. 



3.2. Three-body cluster states in C, B, and C. 

One of the typical examples where cluster and shell features coexist is ^^C. The 
ground state of ^^C is an admixture of Sa-cluster and p3/2-shell closure structures. 
One the other hand, a variety of 3a-cluster states have been suggested in excited 
states in many 3a model calculations since 1970's.^^^ Recently, Tohsaki et al. have 
proposed a new concept of cluster structure in the second 0^ state of ^^C, where three 
a clusters are weakly interacting like a gasMIS Because of the bosonic behavior of 
a particles in a dilute 3a gas state this phenomenon has been discussed in relation 
with Bose-Einstein Condensation of a particles in a dilute nuclear matterHH It is 
a challenging issue to search for such dilute cluster gas-like states in other nuclei, 
for instance, -'^^B and ^^C. Another interesting problem to be solved is a linear-chain 
3a structureM} It has been a long-standing problem whether a linear-chain 3a 
state appears in excited states of ^^C or not. Possibility of linear-chain structures in 
neutron-rich C isotopes is also attracting a great interest as it might be stabilized 
by valence neutrons. 

Applying the AMD method to ^^C, ^^B, and ^^C, we have found various three- 
body cluster structures in their excited states as well as shell-model structures in 
low-lying states. Some are weakly interacting three-cluster states, and some shows 
rather geometric configurations of clusters. We here discuss cluster aspects of these 
nuclei based on the AMD-VAP calculations for ^^C and iiBl2i''iIll'[llF and the /3-7 
constraint AMD calculations for 14qJ59} 
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Fig. 4. Schematic figures for cluster states suggested in ^°Be, ^^Be, and ^^Be. For the molecular- 
orbital states, 2 a cores and the valence neutrons in the a orbital are illustrated. The experi- 
mental values of the excitation energies are also shown. 

3.2.1. Cluster structures of ^'^C 

To investigate excited states of ^^C we performed the AMD-VAP calculations. 
As already mentioned, in the AMD model, all nucleons are independently treated 
without assuming existence of any clusters. The formation or breaking of shell 
structure and clusters is expressed in the twelve-nucleon dynamics after the energy 
variation. Even though any clusters are not a priori assumed, a clusters are formed 
in many excited states of ^^C. Indeed, a variety of 3a cluster states were obtained in 
the AMD results. The experimental and calculated energy levels of ^^C are shown in 
Fig. [6| and density distribution of intrinsic wave functions for the ground and excited 
states is shown in the upper row of Fig. [7| It should be stressed that this is the first 
calculation that succeeded in simultaneously reproducing the energy spectra of the 
ground band and those of excited states having developed 3a cluster structures in 
^^C. The success owes to flexibility of AMD wave functions which can describe both 
shell-model and cluster structures. The calculation also reproduces well transition 
properties such as E2, monopole, and GT transitions. In the result, the ground state 
shows the p3/2-shell closure configuration with a mixing of the Sa-core component, 
while the second 0^ state has a well-developed 3a cluster structure. Since the O2" 
wave function has large overlap with various 3a configurations (amplitudes of the 
superposed wave functions do not concentrate on a specific AMD wave function but 
they fragment largely into various AMD wave functions), this state is regarded as 
the 3a gas-like state, where three a clusters are rather freely moving in a dilute 
density. The 3a chain-like state was suggested in the third 0^ state around 10 MeV 
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Fig. 5. Density distributions of the intrinsic states for the band-head states of ^"^Be and ^^Be 
obtained by AMD-VAPEH'IHI The integrated densities of matter, proton and neutron densities 
are presented in the left, middle and right panels. 



in the AMD result. It is not a "linear" chain but it shows an obtuse triangle 3a 
configuration as seen in the figure. These results are quite similar to those calculated 
with FMD with the unitary correlation operator method (UCOM)© 
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Fig. 6. Left: Energy levels of ^^C calculated with AMD-VAPHl The theoretical levels of ^^C 
calculated with the 3aRGM^^* and 3aGC^^ are also shown. Right: Energy levels of ^^B 
calculated with AMD-VAP.^^' 

3.2.2. Cluster structures of ^^B 

We applied the same method, AMD-VAP, to ^^B and its mirror nucleus ^^C, 
and investigated the structures while focusing on cluster aspect. The energy levels 
of negative-parity states of ^^B were reasonably reproduced by the calculations as 
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Fig. 7. Density distributions of intrinsic states for the ground and excited states of C and B 



calculated by AMD-VAPF*^ 



shown in Fig. [6j 

The 3/23 states at the excitation energy -8 MeV, ^^6(3/23, 8.65 MeV) 
and ^^C(3/2^, 8.10 MeV) are experimentahy known to have abnormal properties 
such as weak GT and Ml transitions compared with normal low- lying states. It is 
also interesting that the 3/2^ state is missing in shell- model calculations for ^^B. 
These facts suggest that the 3/2^ may not be ordinary shell- model-like state but 
may have a developed cluster structure. The AMD calculations for ^^B give good 
results for the energy levels including the 3/2^ state (Fig.[6]) and reproduce well the 
experimental values of transition strengths. For the 3/2^ state, the quenchings of 
GT and Ml transitions are understood because of the developed cluster structure 
of the 3/2^ state which has small overlap with the low- lying shell-model-like states. 
Indeed, the 3/23" states of ^^C and ^^B exhibit the remarkably developed 2a+^He 
and 2a+t clustering (Fig. [7|. 

3.2.3. Analogy of cluster aspects of -'^-'^B to ^'^C 

Comparing the results for ^^B with those for ^^C, we found good analogies of 
cluster aspects between ^^B and ^^C. As shown in Fig. [7| the ground state of ^^B is 
described by the p3/2-shell configuration with a mixing of cluster structure as well as 
that of ^^C. The development of the 2Q:+t-cluster core in the ^^B(3/2^) shows a good 
analogy to that of the 3a-cluster core in the ^^C(2^). The remarkably developed 
2a+t-cluster structures in the ■^■'^B(3/2^) and -^-'^B(l/2^) can be associated with the 
developed 3a cluster in the -^^0(02") and -^^0(03"), respectively. 

Particular attention is paid to analogy of -'^-'^B(3/23 ) to -'^^C(02"). Similarly to the 
case of ■'^^0(02"), the ■'^■'^B(3/2^) wave function has large overlap with various 2a+t 
configurations indicating that the state has no geometric cluster configuration but 
it should be regarded as a weakly interacting 2a+t-cluster state. The root-mean- 
square radius (r.m.s.r.) of the ■'^■'^B(3/2^) state is 3.1 fm and it is remarkably large 
compared with that of the ground state (2.5 fm). Considering amplitudes of the 
wave function fragmented on various configurations and the large radius, the 3/2^ 
state may be a 2Q:+t-cluster state with a dilute density like a gas, where clusters 
are rather freely moving. It should be noted that, for the ^^B(3/2^), the energy 
position relative to the three-cluster break-up threshold is lower and the nuclear size 
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is smaller than those for -'^^C(02"), and therefore the gas-like feature of the ^^B{3/2^) 
might be weaker than the ■'^^C(02"). 

Another analogy between ^^B{3/2^) and ■'^^C(02") is remarkable monopole tran- 
sition strengths from the ground state. The calculated iso-scalar monopole strength 
B{ISO) for the transition 3/2^ ^ 3/2^ is 94 fm^ and the value is in good agreement 
with the experimental value 94zbl6 fm^JH The B{ISO) value is as large as that for 
the monopole transition 0^ O2 in ^^C. 

Recently, structures of -^-'^B has been investigated also by 2a-\-t orthogonality 
condition model (OCM).I^ In the results of 2a+t OCM, different cluster features 
between ^^6(3/23) and ^^0(62") have been pointed out from the point of view of a 
condensation. Cluster gas features of ^^B are under discussion. 

3.2.4. Cluster structures of ■'^^C 

As mentioned, the straight-line chain structure of three a clusters may not be 
stable in ^^C even though the linear- chain- like 3a state with the obtuse triangle 
configuration might exist in the O3" state. We here consider cluster structures of ^"^C, 
in which 3a-core structures with additional two neutrons are expected. 

Structures of excited states of ^^C were investigated with a method of /3-7 con- 
straint AMD in combination with GCM by Suhara and one of the authors (Y. K- 
E.)Il2f We stress again that existence of clusters was not assumed in the model but 
dynamics of fourteen nucleons was solved in the AMD model space. The results 
suggested a variety of developed 3Q;-cluster core structures in excited states. One 
of the new findings is that a 3a linear-chain structure with valence neutrons can be 
stabilized in ^"^C and may construct a = 0^ rotational band above the ^^Be+o^ 
threshold energy (Fig. [s]). As shown in the density distributions of protons and neu- 
trons, the linear-chain state indicates a strongly coupling ^^Be+o: cluster structure, 
where an a cluster is sitting on the head of a deformed ^^Be cluster. It was found 
that additional neutrons play an important role to stabilize the linear-chain config- 
uration. Unfortunately, there is no experimental evidence for the linear-chain state. 
The ^^Be+a decay observationi^B® would be helpful to identify it. 

3.3. Spectroscopy and exotic phenomena in the island of inversion 

Neutron-rich Ne and Mg isotopes around the ^ 20 region called "island of 
inversion" are known to have anomalous properties® These anomalies originate in 
the quenching of the = 20 shell gap in the neutron-rich nuclei and it leads to the 
breakdown of the magic number A^ = 20 and large nuclear deformation. This drastic 
change of nuclear shell structure has been intensively investigated mainly focusing 
on the spectral properties of yrast statesl^'[Z2l'IZ30'[P^' ''^^ Recently thanks to 
the development of the experimental technique, the information of the non-yrast 
stateilll'IZlf'llZl'IZll is rapidly increasing to reveal exotic phenomena peculiar to the 
island of inversion. For example, "the coexistence of spherical and deformed shapes" 
and "the coexistence of normal and intruder configurations" have been discussed 
based on the finding of the second excited 0+ states of ^^M^Z^ and ^^Mg/^^ 

AMD combined with GCM is one of the powerful theoretical approaches to 
investigate the non-yrast states of nuclei in the island of inversion as well as the 
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Fig. 8. Energy levels of ^''C. The theoretical data are those in Ref. 

constraint AMD with GCM by using the Volkov No.2 force with B = H = 0.125 
and the G3RS-type spin-orbit force with ui = —uu = 1600 MeV. Proton and neutron density 
distributions for the linear-chain = 0"*" band are also shown. The dotted lines are the 
experimental and theoretical "'"Be+a threshold energies. The figure is taken from Ref. 67). 



We here introduce some recent AMD studies for shape 
coexistence, one-neutron halo and molecular structure in the island of inversion. 

3.3.1. Many-particle and many-hole states and shape coexistence in ^^Mg 

The neutron orbits and the coexistence of normal and intruder configurations 
in neutron-rich Ne and Mg isotopes are most sensitively probed by the low-lying 
yrast and non-yrast states of odd-mass isotopes. The last neutron's orbit, which 
determines the spin and parity of the ground state, is quite sensitive to the nuclear 
deformation. 

As one of such examples. Fig. [9] shows the calculated energy surface, and pre- 
dicted^i) and observed spectrum of 3iMgll3©® ® ® ® ® 11 ^j^^ ^^^^^^ 

face of ^^Mg (Fig.|9](a)) has the low-lying local minima with different nuclear defor- 
mations. By calculating the occupation probabilities of single-particle orbits, these 
local minima are found to correspond to the Oplh, lp2h, 2p3h and 3p4:h neutron 
configurations in ascending order of deformation. Here the particle-hole configu- 
rations are labeled relative to the = 20 shell closure and nuclear deformation 
becomes larger as numbers of particles in pf shell and holes in sd shell increase. 
To predict the ground state configuration, it is essential to perform a GCM calcu- 
lation. In fact, AMD+GCM calculation (Fig. [o] (b)) shows that the ground state 
is strongly deformed and an almost pure 2p3h configuration in which two neutrons 
are promoted into pf shell across the A^ = 20 shell gap to break down the magic 
number. The calculated magnetic moment of the ground state is — 0.91/XAr, while the 
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Fig. 9. (a) Energy surfaces of the positive- and negative-parity states of ^^Mg as function of 
quadrupole deformation parameter /S calculated by AMD. The AMD results are taken from 
Ref. 81 ). (b) Low- lying spectrum of ^^Mg calculated by AMD+GCM compared with the exper- 
imental assignment suggested in Ref. 83). 



observed value is —O.SSij^n^^^^ Since the spherical Oplh states in which the = 20 
magicity is retained give positive magnetic moment, the breakdown of the = 20 
magic number is confirmed without ambiguity. AMD predicts that the ground state 
is followed by 3/2^, 5/2^ and 7/2^ states with 2p3h configuration to constitute 
the rotational ground band due to the strong deformation. Furthermore, due to the 
quenching of the = 20 shell gap, coexistence of three different configurations at 
small excitation energy is also predicted. The lp2h configuration appears as the 
3/2^ and 7/2^ states at very small excitation energies, the strongly deformed 3p4/i 
configuration constitutes the = 3/2~ rotational band starting from 720 keV, and 
the normal Oplh configuration appears as the 5/2^ state at 1.6 MeV. Most of these 
excited states have been observed^ in good agreement with AMD predictions, by the 
measurements of /3-decayspB®'ll3 one proton or neutron knockout react ion JH*© 
and Coulomb excitation. ^'^ Thus, the coexistence of various neutron configurations 
and deformed states is now established. For the recent discussions on other nuclei, 
readers are directed to Ref. 82). 

3.3.2. Neutron-halo with a deformed core in the island of inversion 

As one of the fascinating phenomena in the island of inversion, we focus on 
the one-neutron halo structure of ^^Ne. Recent experiments at RI Beam Factory in 
RIKEN have revealed the large Coulomb breakup cross section^^' and interaction 
cross sectiorS© of ^^Ne, and the p- or 5-wave neutron-halo structure has been 
suggested from the analysis of the Coulomb breakup Usually neutron-halo 
structure has been discussed based on the "spherical inert core + weakly bound 
neutron" models as done for ^He and ^^Li. However, in this case the core nucleus 
^^Ne is located in the middle of the island of inversion and the assumption of the 
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spherical inert core is inadequate. The last neutron of ^^Ne may be coupled to the 
strongly deformed core with broken magic number. Therefore, the analysis based on 
a full microscopic theory is more suitable and necessary. Since a single-particle wave 
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Fig. 10. (a) Low-lying spectrum of ^^Ne calculated by AMD+RGM. Dashed line shows one neutron 
threshold energy, (b) Proton, neutron, and proton + neutron densities of ^^Ne ground state 
( = 3/2 j") calculated by AMD+RGM. (c) Reaction cross section of Ne isotopes calculated by 
double- folding model using the density distribution obtained by AMD (solid line), AMD+RGM 
(open circle for ^^Ne) compared with the observed The results are taken from 
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function is represented by a Gaussian wave packet, AMD cannot describe properly 
the tail part of halo nucleus which shows the exponential damping. This shortage 
is overcome by combining AMD with resonating group method (AMD+RGM). We 
use the RGM-type wave function for ^^Ne, 



^JM = '^CalA[Xal{r)Yim{r)(l)n(l)30^^{a)] 



JM ' 



(3-1) 



al 



where (^n and XQ;K^)^^m(^) are the spin wave function of the valence neutron and the 
relative motion between the valence neutron and ^^Ne, respectively. 030^^(0:) is the 
internal wave function of ^^Ne solved by AMD+GCMp^ and a labels the ground 
and excited states of ^^Ne. Since the AMD wave function of ^^Ne is a superposition of 
Slater determinant of Gaussians, the coefficients c^i and the relative wave function 
Xnlir) can be solved by transforming Eq. (3-1) into a GCM type wave function. 



Recent development of high-performance computing enabled this CPU demanding 
calculation. 



Figure 10 (a) shows the predicted spectrum of ^^Ne by AMD+RGM. The ground 
state is the 3/2~ state and the calculated one neutron separation energy is Sn — 0.45 
MeV, while the observed value is 0.29=bl.64 MeV.^^^ Due to this small separation 
energy, the density distribution of ^^Ne calculated by AMD+RGM shows the long 
tail at large distance indicating a p-wave halo (lp3/2-neutron coupled to ^^Ne) as 
shown in Fig. [To] (b). The effect of the deformed core appears as the core excitation. 
From the coefficients Cc^/, it is found that the p^i2 0^^Ne(2+) configuration amounts 
to 41%, that is larger than the ^3/2 ® ^^Ne(O^) configuration which amounts to 
37%. These values suggest that the P3/2 neutron is coupled to the deformed and 
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rotating ground band of ^^Ne. Using AMD and AMD+RGM wave functions, the 
reaction cross sections of ^^Ne and other Ne isotopes are analyzed based on the 
double-folding Figure 



10 



(c) compares the calculated and observed 
reaction cross section. We can see that AMD wave function shows overall agreement 
with the observation except for ^^Ne, and anomalous large cross section of ^^Ne 
is reasonably described by employing the AMD+RGM wave function. Thus, with 
a help of high-performance computing, AMD combined with RGM and reaction 
theory is a promising method to investigate neutron-halo nuclei and their reactions 
in heavier mass region. We can find many other candidates of weakly bound system 
such as ^^Mg and ^^Mg from the systematics of the binding energy, and analysis of 
them is now ongoing. 

3.3.3. Molecule-like states in the island of inversion 

Another example of exotic phenomena in the island of inversion is the molecule- 
like structure at highly excited region. As discussed in the subsection 3.1, several 
excited states of O, F and Ne isotopes have been predicted to have molecule-like 
structure analogous to Be isotopesll2l'l2§f Especi ally , the candidates of molecule-like 
states have been recently observed in 18~20 q|99} ,[ioo| ,[ioi| found to qualitatively 
agree with the AMD predictions Ml Thus the exploration of molecule-like structure 
is expanding to heavier neutron-rich systems. In the case of F and Ne isotopes, 
since the neutron drip line is farther extended than for O isotopes, we can expect 
molecule-like structure with more valence neutrons and more exotic phenomena. 



14 



„12 

> 

1.10 
^ 8 



a 6 

CD 
-t— ' 

4 



proton hole 
with two neutrons in a2 

1/2- 



3/2+ \l/2- 



proton hole 3/2"*" 
.with one neutron in ct 



proton hole 
with no neutron in c1'2" 
1/2- 



5/2+ 5/2+ 



3/2+ 

1/2^ 
1/2- 



5/2+ 



1/2- 

3/2+ 
1/2^ 



5/2+ 



21 1 



23r 



25r 



27r 



1/2- 

3/2+ 
1/2- 

"5/2+ 



proton hole with ct^ 



proton hole with <j 



29c 






ground k2 




Fig. 11. Left: Band-head energies of proton hole states of F isotopes. Red Unes show the states with 
a proton excitation from p to sd shell and no neutron in pf shell, while green (blue) lines show 
those with one (two) neutron(s) in pf shell that have molecule-like structure. Right: Density 
distribution of the core (^^F, contour lines) and two valence neutrons (color plot) of ^^F. From 
top to bottom, each panel corresponds to the states with two, one and no neutrons in pf shell 
with a proton hole and the ground state, respectively. This figure is taken from Ref. |7|. 



Figure 11 shows the energies of the F isotopes obtained by AMD+GCM^^^ as 
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an example. All of these F isotopes have molecule-like states in their excited states 
and there are always two kinds of molecular bands (green and blue lines in Fig. [IT] 
left). In these states, one or two valence neutrons are excited into a orbital {pf shell) 
together with a proton excitation from p to sd shell, and the a+N cluster structure 



develops (Fig. 11 right). It is notable that the drastic reduction of their excitation 
energies toward the neutron drip line is predicted. To understand this reduction, 
readers are reminded following points. (1) As mentioned in section 3.1, the a orbital 
originates in the pf shell, and its energy is lowered in the island of inversion due 
to the quenching of the = 20 shell gap. (2) If the core has cluster structure, it 
induces the deformation of the system and it further lowers the energy of a orbital 
as in the case of Be isotopes. Therefore, the neutron excitation into a orbital and 
a clustering of the core work in a cooperative way to reduce the excitation energies 
of molecular states in the island of inversion. Up to now, several candidates for the 
molecular states in lighter F isotopes are experimentally known il2H} '1122} '[l2if More 
data for F isotopes near the drip line will be experimentally available in near future. 

3.4. Superdeformation in sd-shell nuclei 

A recent development in high-spin physics is the di scov ery o f the superdeformed 
band at very smah mass region {A - 40) such as ^^Ai^ 40QgfTo6^ "^"^Ti^In 
contrast with heavier nuclei, these superdeformed bands are assigned from high-spin 
states down to very low-spin states (J^ = 0^ ^ 4+). From this discovery, the 
low-lying 0^ states known for a long time are now identified as the band-head of 
the superdeformed bands. On the other hand, a clustering of low-lying states and 
molecular structure of highly excited states have long been discussed. AMD studies 
have been made to investigate properties of superdeformed states and to reveal the 
relationship between the superdeformed states and cluster states. 

3.4.1. Dual nature of superdeformed states and its evolution to the molecular states 
The strongly deformed excited state of ^^S is of particular interest because of 
the following reasons. 1) ^^S {N — Z — 16) is a double mafflc nucleus of superdefor- 
mation. A couple of mean- field calculation^l22} '11121 have predicted the superde- 
formed state around ^ 10 MeV and it has \hbj excited configuration relative to 
the ground state. 2) Using unique optical potential for ^^O-^^O scattering jllll'Ilill 
^^O+^^O cluster model^^^' showed the presence of three molecular bands. Among 
them the lowest energy band is located a few MeV below the ^^O+^^O threshold 
energy and coincides with the predicted superdeformed band mentioned above, while 
the highest energy band nicely reproduces the well-known ^^O+^^O molecular res- 
onances .Hill* 'Eill* Thus the two different theoretical approaches give qualitatively 
the same result and suggest the relationship between the superdeformed states and 
cluster states. 

AMD gives a unified understanding for the results of the mean-field and clus- 
ter model studies and shows that the superdeformed band evolves to the ^^O+^^O 
molecular bands as the excitation energy increases. Similar to the mean-field cal- 
culations, AMD predicts a well pronounced superdeformed minimum (Fig. [12] (a)). 
The density distribution of the superdeformed wave function clearly shows ^^+^^0 
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Fig. 12. (a) Energy surface of ^^S as function of quadrupole deformation /3. Solid lines show the 
energy surfaces for each angular momentum obtained by AMD, while the dashed line shows 
that obtained by ^^O+^^O cluster wave function, (b) Three rotational bands built around the 
superdeformed minimum obtained by the AMD+GCM (sohd lines) and ^^O+^^O cluster model 
wave functions. The principal quantum number N for each rotational band is also shown. The 
figures are taken from Ref. 116). 



molecule-like structure and the distance between ^0+ O increases as deforma- 
tion becomes larger. Note that the superdeformed wave function is smoothly con- 
nected to the ^^O+^^O cluster wave function around the Coulomb barrier (/3 ^1.1) 
in AMD. Around the superdeformed minimum, the ^^O clusters are distorted by the 
formation of the mean-field and by the spin-orbit interaction to gain more binding 
energy, while the pure ^^O+^^O cluster structure is restored as inter-cluster dis- 
tance (deformation) becomes larger. The AMD+GCM calculation has shown that 
three rotational bands appear as the superposition of the wave functions around 
this superdeformed minimum (Fig. [12] (b)). The lowest energy band coincides with 
the superdeformed band predicted by the mean- field calculations. The ^^O+^^O 
cluster component in this band amounts to 57% which means the superdeformed 
band has the ^^O+^^O molecule-like structure, but ^^O clusters are considerably 
distorted. This shows dual nature of this band (mixing of mean-field and cluster 
nature) and explains why the mean-field calculations and cluster model give quali- 
tatively the same results. The degree-of-freedom of inter-cluster motion embedded 
in the superdeformed band shows up as two excited rotational bands. The cluster 
components of these two bands amount to more than 90% and these bands are in- 
terpreted as the excitation mode of the superdeformed band in which the relative 
motion between two ^^O clusters is excited by 2 and 4 huj^ that are confirmed by the 
analysis of the inter-cluster motion of the AMD wave function. The highest band 
with 4:h(jj excitation of relative motion plausibly agrees with the observed ^^O+^^O 
molecular bandlHlf'lH^ Thus the superdeformed band and the molecular band can 
be regarded as a series of the -'^^O+^^O cluster bands, and the superdeformed band 
evolves to the ^^O+^^O molecular band as the inter-cluster motion is excited. 

3.4.2. Superdeformation and clusters in A ^ 40 region 

The superdeformed state and ^^O+^^O clustering of ^^S shed light on the rela- 
tionship between the (super) deformed states and clustering. As an illustrative and 
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interesting example, Fig. [Ts] (a) shows the observed spectra of doubly magic nucleus 
^^Ca. Two deformed rotational bands start from the 0+ states at 3.35 and 5.21 
MeV. It has been suggested that the former band is predominated by a Ahuj excited 
configuration and the latter is a Shuj configurationlliZi '1112} The discovery of the 
high-spin states^^^^ identified the latter band as a superdeformed band, and several 
theoretical studies based on the mean-field models have been performedli2§'lii2l On 
the other hand, based on the clust er m odel, a+^^Ar clustering of the band starting 



from 3.35 MeV has been suggestedli2f'llH2^ and experimentally confirmed by the a 
transfer reaction on 36^^ I2i),i22l ^^le following, we denote the band starting at 
3.35 MeV as ND band and the band at 5.21 MeV as SD band. 
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Fig. 13. (a) Observed positive-parity states of ^°Ca. (b) and (c) Calculated positive- and negative- 
parity spectra by AMD+GCM.^^^' 'OB', 'ND', 'SD' and 'N3' denote the oblate and prolate 
deformed bands, superdeformed band and parity doublet partner of the superdeformed band, 
respectively. (dl)-(d9) Intrinsic density distribution obtained by the constraint variational cal- 
culation (see text), (e) Calculated spectra of the highly excited positive-parity states. HNl 
and HN2 denote the higher nodal ^^C+^^Si cluster bands in which the inter-cluster motion is 



excited. The AMD results are taken from Ref. 123) 
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AMD studies have been performed to clarify nature of the ND and SD bande^H^'llHlf 
to reveal the coexistence of the different deformed bands and their relationship to 



the cluster structure. Here we mainly discuss the result reported in Ref. 123), in 
which the variational calculation is performed under the constraint on the principal 
quantum number of the harmonic oscillator, 



i=i i=i 



(3-2) 



Here the lowest Pauli allowed value is = 60, and the average number of nucleons 
promoted from sd to pf shell increases for a larger value of N. Figure [Ts] (dl)-(d9) 
shows the intrinsic density distributions obtained for different values of N. For a 
small value of N, an almost spherical state is obtained (dl), and prolate (d2) and 
oblate (d3) deformed states appear by slight increase of N. Further increase of N 
develops the parity asymmetric prolate deformed state (d4). The parity asymmetry 
of this wave function originates in ^^C+^^Si cluster nature of the superdeformed 
state, since it evolves into a prominent ^^C+^^Si cluster state as N increases (d4)- 
(d6). An important point is that similar -'^^C+^^Si cluster-like state (d9) also appears 
in the negative-parity state with large A^, while the negative-parity state with small 
A^ (d7) corresponds to the Ihw excited state built on the spherical state (dl). 

Several rotational bands are obtained as the superposition of those wave func- 



tions by AMD+GCM as shown in Fig. 13 (b) and (c) compared with the observed 



data (Fig. [13] (a)). In the positive-parity states, there are two rotational bands (de- 
noted as ND and OB in Fig. [13] (b)) that are dominantly composed of the wave 
functions (d2) and (d3) respectively, predicting a prolate and oblate shape coexis- 
tence. Though the calculated excitation energy is larger than the observation, the 
calculated ND band is assigned to the observed band starting from 3.35 MeV from 
the comparison of their B{E2) strengths, while the corresponding oblate deformed 
band has not been experimentally assigned yet. 

The highlight of the result is the superdeformed band and related ones that are 



associated with the parity-asymmetric -'^^C+^^Si clustering nature (Fig. 13 (d4)-(d6) 
and (d8)-(d9)). The calculated superdeformed band (SD in Fig. [13] (b)) is dominated 
by the wave functions (d4)-(d6) and starts around 8 MeV, while the observed band is 
located at 5.21 MeV. The AMD calculation has predicted the presence of the parity 
doublet partner of the superdeformed band generated by the parity asymmetry. The 



— band starting around 11 MeV in the negative parity (N3 in Fig. 13 (c)) 



is mainly composed of the wave function shown in Fig. 13 (d9) that also manifests 
i2(-i_^28gj clustering and can be understood as the parity doublet partner of the 
superdeformed band. The clustering nature of the superdeformed band also appears 
as the two excited molecular bands that are predicted at high excitation energy 
region (HNl and HN2 in Fig. [Ts] (e)). These bands are dominated by the wave 
functions with A^ larger than for the superdeformed band dominated by (d4)-(d6), 
which indicates that they are the nodal excitation modes of the inter-cluster motion. 
Similar to the case of ^^S, the superdeformed band (SD) and nodal excited bands 
(HNl and HN2) can be understood as a series of the ■'^^C+^^Si cluster bands. Thus, 
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Fig. 14. Schematic figure summarizing the excitation modes of ^°Ca. Prolate and oblate deformed 
states coexist at small excitation energy. The superdeformed band has dual nature of the mean- 
field and ^^C+^^Si clustering. The triaxiality of the superdeformed band generates the = 2^ 
sideband, while parity asymmetry of the ^^C +^^Si cluster structure generates the parity doublet 
partner. Furthermore, the nodal excitation of the inter-cluster motion generates a series of the 
i2q _|_ 28gj molecular bands. 



the observation of the parity doublet partner and nodal excited bands will be essential 
to reveal clustering nature of the superdeformed band. It is also noted that the 



triaxial deformation of the superdeformed band has been also discussed in Ref. 124), 
and it has been suggested that the = 2+ band accompanies the superdeformed 
band. 



The situation of the ^^Ca is schematically summarized as shown Fig. 14, Above 
the ground state, a prolate and oblate shape coexistence presents as two rotational 
bands. The superdeformed band is also located at the similar energy region and 
manifests the dual nature of the strongly deformed mean-field and ^^C+^^Si clus- 
tering. The AMD calculation has predicted that the ^^C+^^Si clustering nature of 
the superdeformed band generates the parity doublet partner in the negative-parity 
states from the parity asymmetry and the nodal excited bands at highly excited 
region from the excitation of the inter-cluster motion. 

3.5. Decoupling between Proton and neutron deformations 

For Z N unstable nuclei, exotic phenomena concerning nuclear deformation 
have been discovered. If the shell effect for proton orbits and that for neutron 
ones compete with each other, the shape of proton density may be affected by the 
neutron structure, or it might be insensitive to neutron structure. In the latter case, 
decoupling of deformations between proton density and neutron one may occur. 
The decoupling, i.e., different proton and neutron deformations is possible in light- 
mass nuclei and it can be observed in the quadrupole transition properties such as 
the ratio of the neutron transition matrix amplitude to the proton one (so-called 
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Mn/Mp ratio). 

Such a decoupling between proton and neutron shapes has been suggested, for 
instance, in ^^C for which an enhanced Mn/Mp ratio cahed neutron dominance has 
been observed in the ground-band transition, 2^ 0^, in inelastic scattering. 
The neutron dominance was described by opposite deformations, i.e., a n oblate pro- 
ton shape and a prolate neutron one (Fig. 15) in the study with AMD1125}'112ZI The 
opposite deformations in ^^C has been supported also by the abnormally small E2 
transition strength B{E2;2^ 0^) =2.6 ±0.9 e^fm^ observed by the life time 
measurement of the i^C(2+)Ill§ The strength S(£;2;2+ 0^) in ^^C is smah 
compared with those for other C isotopes, ^^C, ^^C, and ^^C. The hindrance of the 
B{E2; 2^ 0^) implies a small proton deformation, however, it seems to contradict 
a large deformation expected from the small excitation energy of the 2^ state if 
proton and neutron deformations are assumed to be consistent as usual. 

We studied the structure of even-even C isotopes with AMDll2£f'll22f The sys- 
tematics of the binding energies, radii, and E2 transition strengths of C isotopes 
were qualitatively reproduced by a simple version of AMD calculations (VBP). The 
results indicate that the neutron shape drastically changes depending on the neu- 
tron number, while the proton shape is rather stable and insensitive to the neutron 
structure. One of the striking features is that the difference between proton and 
neutron shapes was suggested in ^^C and ^^C in which prolate neutron shapes are 
favored. In spite of the prolate neutron structure, the proton structure shows an 
oblate deformation resulting in the opposite deformations. The deformation param- 
eters for proton and neutron densities of the intrinsic state are 7^) = (0.41, 0.277r) 
and (/3n,7n) = (0.53, O.OOtt) for ^^C, and they are (/3p,7p) = (0.32, 0.267r) and 



(/3n,7n) = (0.34, O.OOtt) for ^^C (Figs. 15 and 16). The reason for opposite pro- 



ton and neutron deformations is that a Z = 6 system favors an oblate proton shape 
because of the proton shell effect while a = 10 or = 4 nucleus has prolate 
trends of the neutron shape due to the neutron shell effect. In other words, the 
Z = 6 proton structure is not so much affected by the neutron structure but it keeps 
the oblate tendency. 

To discuss the neutron deformation, mirror analysis is useful. In the mirror 
analysis for ^^C and ^^Be, the neutron transition matrix Mn for the ground-band 
transition is evaluated from B{E2) in ^^Be by assuming mirror symmetry. The 
experimental value of the M^/Mp ratio in ^^C deduced by the mirror analysis is 
described by the AMD calculation, and it can be understood with the opposite de- 
formations between proton and neutron densities (Fig. [16]). The neutron dominance 
in the ground-band transition is more remarkable in as seen in the theoretical 
results. Unfortunately, there is no direct data of the E2 strength for the mirror 
nuclei of ^^C, however, as mentioned before, the observed inelastic scattering cross 
section implies an enhanced Mn/Mp ratio indicating the neutron dominance ilH5f It is 
worth to mention that microscopic coupled-channel calculations with the transition 
densities obtained by the AMD calculation have reproduced the inelastic scattering 
data successfullylHZf 

To clarify the oblate shape of the proton structure in ^^C and ^^C, observations 
of possible K — 2 side bands and their transition properties would be helpful probes. 
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Such a nucleus having oblate proton and prolate neutron structures may show an 
isovector triaxiality. If the case, a K — 2 side band can be constructed from the 
rotation around the symmetric axis of the prolate neutron part instead of the rotation 
around the perpendicular axis for the K — {) ground band (Fig. 16). Since the proton 
contribution should be dominant while the neutron contribution is minor in the 
rotation for the K — 2 side band, the inter-band transition 22 0^ may show the 
proton dominance resulting in a small M^/Mp ratio. In fact, the calculated Mn/Mp 
ratios for the 2^ 0^ in ^^C and ^^C are quenched and they indicate the proton 
dominance. There is no experimental information of the transition strength for 
22" ^ 0^. Inelastic scatterings of ^^C and ^^C will be good probes to experimentally 
confirm the proton dominance. 



(a) i^C (c) i^C 

, neutron ^ J.=K=0 
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Fig. 15. Schematic figures for different shapes of proton and neutron densities in (a) ^°C, (b) ^°Be, 
and (c) ^^C. (d) Surface cut at constant proton and neutron densities of ^^C obtained by the 
VBP calculations with AMD. (e) Prolate neutron density of ^^C. (f) Oblate proton density of 



§4. Applications of time-dependent AMD method to nuclear response 

and reaction 

Interesting phenomena, including clustering phenomena, appear in nuclear many- 
body systems not only near the ground state but also in excited states in wide range of 
excitation energies. These include giant resonances of collective modes and stronger 
expansion of the system which is often followed by the disintegration into fragment 
nuclei. Not only the collective degrees of freedom but also non-collective (thermal) 
excitations play important roles in many cases. Since the number of states increases 
very rapidly as the excitation energy goes up, it is difficult to describe these phe- 
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Fig. 16. Left: the experimental energy levels of the 0^, 2^, and 2 J states and Mn/Mp ratios for 
the 2+ of and 2+ of and ^^C^\^ The experimental values of the neutron 

matrix amplitude (Mn) are deduced from the corresponding B(E2) values of the mirror nucleus. 
The values in the parentheses are the theoretical values for Mn/Mp of the AMD calculation UHl 
Right: deformation parameters for the intrinsic wave functions of and ^^C. The filled trian- 
gles indicate /Sp and jp for the proton part and the open circles are /3n and jn for the neutron 
part. 



nomena based on individual quantum states. The time-dependent version of AMD 
has been a powerful approach for such problems of highly excited systems. 

4.1. Dipole resonances 

Among exotic phenomena in neutron-rich nuclei, one of the hot subjects is dipole 
strengths in low energy regions, which are expected to enhance in neutron-rich nuclei 
due to excess neutrons. For instance, soft resonances in extremely low-energy regions 
and pygmy resonances below giant dipole resonance (GDR) energy are attracting 
great interests (see, for example, Refs. |132[ ), 133), 134), |135 ), |136 ), 137| ) and references 
therein). To investigate isovector dipole responses of neutro n-ric h nuclei, we applied 
a time-dependent version of AMD without stochastic termslHH 

In order to calculate responses to external fields, we first solve the static problem 
and obtain the optimum parameter set Z^, which gives the energy minimum state 
^amd(^°) in the AMD model space. Then, we boost the ^amd(^^) instantaneously 
at t = by imposing an external perturbative dipole field. 



^(t = 0+) = e-^^^#AMD(^°), (4-1) 

VMr,t) = eF{r)5{t), (4-2) 
A 

F{r) = M{El,fi) = J^e-V,yi^(fi), (4-3) 



where e is an arbitrary small number and e^^^ is the El recoil charge, Ne/A for pro- 
tons and —Ze/A for neutrons. Note that the initial state ^{t = 0+) after imposing 
the dipole field is written with a single AMD wave function ^amd(^(^ = 0+))- 

Following the time-dependent AMD method, we can calculate the time evolution 
of the system, ^{t) = ^AMD{Z{t)), from the initial state ^{t = 0+) by using the 
equation of motion Eq. ( |2-5[ ). Once the wave function ^{t) is obtained as a function of 
time, the transition strength can be obtained by Fourier transform of the expectation 
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value of M.{E1, fi) as follows, 

- Y.\{n\M{El,m?5{u: - u^) (4-4) 



duj 

n 

= Im / dt{^{t)\M{El,fi)\^{t))e'^\ (4-5) 

Jo 

where |0) is the ground state and \n) is the excited state with the excitation energy 

In the present framework, dB{(jj; El) /dcj consists of discrete peaks in principle, 
because the present AMD method is a bound state approximation and continuum 
states are not taken into account. We introduce a smoothing parameter F by hand 



in the Fourier transform in Eq. (4-4) which may simulate the escape and spreading 
widths of resonances. 

One of the advantages of time- dependent AMD is that we can obtain intuitive 
interpretations for each mode by analyzing time evolution of Gaussian centers Z^(t) 
of single-particle wave packets. Another advantage is that the present method is 
free from the spurious center-of-mass motion because center-of-mass motion can be 
exactly separated from ^{t). 

We applied this method to Be, B, and C isotopes and investigated the El les- 
onances.^^^^ The El strengths are shown in Fig. [l7| It was found that remarkable 
peaks appear in ^^Be, ^^B, and ^^C in the E* = 10 — 15 MeV region decoupling 
from the GDR. Those soft dipole resonances arise from the relative motion of ex- 
cess neutrons against a core, which is decoupled from the motion inside the core. 
In other words, the soft resonances appear due to the excitation of excess neutrons 
around the rather hard core. In fact, the strengths of the soft dipole resonances 
almost exhaust the cluster sum rule values for the core and valence neutrons. In 
further neutron-rich B and C isotopes with N > 10, the strengths for the soft dipole 
resonances decline compared with those in ^^B and ^^C. The reason for decreasing 
low-lying strengths is that motion of excess neutrons assimilates into neutron motion 
inside the core and its decoupling from the core weakens. As a result, the excita- 
tion energies of the GDR decrease with the enhancement of the neutron skin. It is 
striking that strengths of the soft dipole resonances do not necessarily increase with 
the increase of excess neutrons. Instead, the feature of the soft resonances rapidly 
changes depending on the proton and neutron numbers of the system. The strengths 
of the soft dipole resonances depends on how much the coherent motion of the excess 
neutrons decouples from the motion inside the core. 

4.2. Radial oscillations 

Radial breathing oscillations of nuclei are interesting particularly because they 
are closely related to the compressibility of nuclei and nuclear matter. The AMD 
approach is, however, not very suitable for the precise analysis of the monopole 
strengths since the state boosted by the monopole operator from the ground state 
AMD wave function is no longer an AMD state. Nevertheless, AMD can be a power- 
ful tool to explore the radial oscillations at various amplitudes. The large amplitude 
oscillations should be continuously linked to the expansion of nuclei without restora- 
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Fig. 17. El transition strengths (e^fm^/MeV) of Be, B, and C isotopes calculated with the time- 
dependent AMD method by using the MVl (m = 0.576 and 6 = /i = 0) + GSRS {m = -un = 
900 MeV) force. ^"^^^ The smoothing parameter is chosen to be = 2 MeV. Thin dash-dotted, 
solid, and dotted lines are the contribution of vibration for the x,y, and z-directions, respectively. 
The total strengths are shown by the thick solid lines. 



tion at higher energies where many-body correlations play important roles in forming 
clusters and fragment nuclei. Therefore the possibility of cluster correlations in lower 
energy oscillations is an interesting question. 



20 
15 
10 

^30 

CD 

1.25 

3 

30 
25 
20 





SLy4 - 
Sill -- 
-Gogny ■■■ 





















10 15 
E* [MeVl 



20 



Fig. 18. AMD prediction for the frequency huu of the radial oscillation as a function of the oscillation 
energy for ^^C, ^^O and ^°Ca with the SLy4, SIII and Gogny interactions. The nuclear matter 
incompressibilities for these interactions are K = 230, 355 and 228 MeV, respectively. The 



figure is taken from Ref. 139). 



The radial oscillations of the ^^C and other nuclei were calculated with AMD [Eq. 
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(2-5)] without any stochastic terms in Ref. 139). The time evolution is solved from 
the initial state prepared by placing three a clusters on a regular triangle in the ^^C 
case. By changing the size of the initial regular triangle, oscillations with different 



amplitudes were studied. Figure [T8| shows the dependence of the oscillation frequency 
on the amplitude or the excitation energy. The results for three different effective 
interactions are shown. Dependence on the incompressibility is clearly observed. 
The dependence on the amplitude suggests the unharmomicity of the oscillation. It 
should be noted that the one-phonon excitation corresponds to the excitation energy 
oi E"" — hjjj where uj is the angular frequency. 

The calculations of Fig. [18] were carried out with the width parameters v that op- 
timize the ground state energy. Depending on the chosen width parameter, however, 
the frequ ency sometimes shows anomalous behavior as in the system of a double- well 
potentials 
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Fig. 19. (a) Fourier transform FT[r](c(;) of the FMD oscillation pattern in small amplitude, with 
the same initial condition as in the AMD calculation, (b) The peak position /icj as a function 
of the oscillation amplitude. The figure is taken from Ref. |139|). 



Th e sa me problem was studied by the fermionic molecular dynamics (FMD) 
as well 11221 In FMD, the width parameters z^^(t) of individual wave packets are 
treated as time-dependent variables as well as the centroid variables. In this case, the 
calculated result shows an oscillation pattern composed of two modes with different 



frequencies as shown in Fig. [19) By analyzing the motions of the variables, it was 
found that one of the two modes corresponds to the change of the wave packet 
widths and the other corresponds to the motion of the wave packet centroids. The 
latter mode is the same as that observed in the AMD result and related to the 
^-clustering degrees of freedom. The former mode is the breathing of individual 
single-particle wave packets. Thus these calculations suggest that the single-particle 
excitation and the clustering excitation are both important in radially oscillating 
systems. The simplest version of AMD without any stochastic terms can describe 
clustering excitations. 

4.3. Multifragmentation in expanding systems 

It has been well known experimentally that a lot of fragment nuclei are pro- 
duced in each event of heavy-ion collisions in various situations if the incident en- 
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ergy is more than ten MeV/nucleon. Multifragmentation is an interesting problem 
in excited nuclear systems, in which quantum many-body correlations play essential 
roles as well as the existence of nuclear liquid-gas phase transition. The microscopic 
description of multifragmentation is, in principle, a highly complicated problem of 
quantum many-body systems. Transport models have been developed for heavy-ion 
collisions with some classical approximations. Compared with other transport mod- 
els, some quantum features have been incorporated into AMD by employing fully 
antisymmetrized wave functions. 




Fig. 20. Time evolution of density profiles obtained by the AMD (upper) and SMF (lower) models, 
for the central collision of ^^^Sn + ^^^Sn at 50 MeV/nucleon. 



Figure 20 shows a typical example of the time evolution of ^^^Sn + ^^^Sn central 
colhsions at 50 MeV/nucleonIll2l'El} The calculation by AMD in the upper row 
is compared with the calculation by the stochastic mean-field (SMF) model in the 
lower row. In both cases, the system compressed at an early stage starts to expand 
almost spherically. As the expansion proceeds, the density fluctuation develops to 
form many fragments. 

The AMD calculation was performed with the coherence time r = tnn- The 
SMF model, which is based on the single-particle motion in the mean field, takes 
into account the two-nucleon collisions and fluctuations*112^ Therefore these two 
models are conceptually similar, but the results can be different due to the different 
approximate treatments of fluctuations. In fact, it is observed that the density 
fluctuation (among different events) is already developing in AMD at the relatively 
early stage of 50 < t < 100 fm/c, while the fluctuation develops in SMF only 
at a later staget ^ 100 fm/c suggesting a fragmentation mechanism by spinodal 
decomposit ionium Thus the many-body correlations are stronger in AMD. This 
difference can be interpreted as the origin of the differences in the e xpansion velocity, 
the nucleon emission and so on predicted by these models Ill2l 'Hill 

The final results of these models for the fragment charge distribution have been 



compared with the experimental data with reasonable successes llllf'lHHf Figure 21 
shows the AMD results compared with data. The result depends very much on the 
choice of the coherence time. In this reaction system, the fragment yields for Z ^ 3 
are well reproduced by AMD when the wave packet splitting with the coherence time 
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Fig. 21. The charge distribution of the produced clusters in ^^^Xe+Sn colhsions at 50 MeV/nucleon 
with the impact parameter < b < 4 fm, after calculating the secondary decay of excited clusters 
and applying the experimental filter for the detector setup. Solid histogram (labeled AMD/DS) 
shows the result of AMD with the coherence time r = tnn, while the dotted histogram (labeled 
AMD/D) shows the result with the strongest decoherence r ^ 0. The INDRA experimental 
data are shown by solid points. The figure is taken from Ref. 29). 



^ — '^NN is introduced. It is often convenient to define the liquid and gas parts of the 
system as the parts composed of Z > 3 fragments and Z < 2 particles, respectively. 
The comparison shows that the total charge of the liquid part, Ziiq, in the AMD 
result (with r = tnn) is consistent with the experimental data, and therefore the 
total charge of the gas part, Zgas = ^system — ^liq, is also consistent. 

However, a problem is found in the composition of the gas part in the result with 
^ = "^NN- The a-particle multiplicity ^ 7 is too small and the proton multiplicity 
Mp ^ 20 is too large compared with the experimental data Ma ^ Mp ^ 10. It should 
be noted that only about 10% of the total protons in the system is emitted as free 
protons in this reaction at 50 MeV/nucleon. It is also known experimentally that 
still a half of the protons are bound in clusters even at 1 GeV/nucleonEll Thus 
the experimental data have been suggesting the importance of cluster correlations. 
The comparison with the data indicates that the AMD approach with the wave 
packet splitting, which is largely based on the single-particle motion in the mean 
field, does not include sufficient cluster correlations in the dynamics. The problem 
of the gas composition may influence on the liquid part of the system because the 
energy balance and the number of the effective degrees of freedom will change as the 
gas composition changes. Therefore the proper treatment of cluster correlations in 
dynamical approaches is an urgent issue. 
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Fig. 22. Isotope distribution in ^^C + ^^C reaction at 28.7 MeV/nucleon. Circles are the yields 
before the secondary decay, while triangles are those after the secondary decay. Squares are the 
experimental data. Lines connect isotopes. The figure is taken from Ref. [2|. 



4.4. Fragmentation in collisions of light nuclei 

In an early study with AMD, the fragmentation in the ^^C + ^^C reaction was 
studied at 28.7 MeV/nucleonllf'l2f The fragment isotope distribution is reproduced 
weh by AMD as shown in Fig. [22j Especiahy the large production cross section of a 
particles is well reproduced in this case. The calculation was done with the stochastic 
two-nucleon collisions but without wave packet splitting. Thus the situation here is 
different from the multifragmentation in heavier systems as seen in the previous 
subsection where the wave packet splitting is very important. 

The detailed mechanism of the fragmentation of the ^^C projectile was studied 



in Ref. '146) by Takemoto et al. The solid line of the lower part of Fig. 23 shows the 
distribution of the excitation energy of the ^^C projectile after the interaction with 
the target in ^^C + ^^N reaction at 35 MeV/nucleon. In addition to the big bump 
around ^ 30 MeV, there is a sharp peak at ^ 10 MeV. The a-clustering 
states of ^^C in this energy region are strongly excited by the heavy-ion reaction, 
and they contribute to the large yield of a particles. In fact, the dashed line shows 
the contribution from the events where ^^C broke up into three a particles. On the 
other hand, in the proton induced reactions, the excitation of ^^C is mainly of single 
particle nature and the a clustering states are not excited at all as shown in the 



upper part of Fig. 23 



For the multifragmentation in ^^Ca + ^^Ca collisions at 35 MeV/nucleon, the 
wave packet splitting plays important roles to enable the breakup of the system into 
small pieces. Figure [24] compares the results without wave packet splitting and with 
the wave packet splitting in the limit of the small coherence time r ^ (i.e., the 
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Fig. 23. The distribution of excitation energy of ^^C projectile before its breakup. Upper figure is 
for p + ^^C collisions, while the lower figure is for ^^C + ^^N collision at 35 MeV/nucleon. The 
figure is taken from Ref. 146). 



strongest decoherence) . The latter result is consistent with the experimental data. 
Without wave packet splitting, the two nuclei go through each other without forming 
fragments of intermediate size. The wave packet splitting allows the mixing of the 
two nuclei or the neck formation, so that more than two fragments can be formed 
from the system expanding in the beam direction. The a-particle multiplicity also 
depends on the wave packet splitting very much. 

In short, the strength of the wave packet splitting (viewed from AMD) or the 
decoherence of the single-particle states (viewed from the mean-field theory) is a key 
ingredient for the description of fragmentation. For the stronger splitting, the system 
tends to expand strongly and to break into small fragments and many a particles. 
Unfortunately, the appropriate strength seems to depend on the size of the system 
and/or the incident energy. A more consistent understanding may be possible if the 
cluster correlations in dynamical systems are more explicitly treated. 

4.5. Statistical properties of excited systems 

One of the aims of the study of multifragmentation has been to extract the 
information of excited nuclear matter in which the liquid-gas phase transition is 
expected as in the system of the Van der Waals equation of state. The concept of 
phase transition in finite many-body systems, a s in the heavy-ion collision systems, 
has been improved very much in recent studies.^^^'EH* Phase transition is clearly 
defined in finite systems by considering microcanonical ensembles. 

By solving the time evolution of a many-body system in a container of a finite 
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Fig. 24. Mass distribution in '^°Ca+'^°Ca collision at 35 MeV/nucleon. Calculated results of AMD 
before and after the secondary decay are shown by dotted and solid histograms, respectively. 
The left figure is the result by AMD without wave packet splitting, and the right figure is the 
result with the strongest decoherence r ^ 0. The figure is taken from Ref. 28). 



volume for a very long time, it is possible in principle to generate a microcanonical 
ensemble for the given energy and volume. However, it is a non- trivial question 
whether dynamical models such as AMD can produce a correct statistical ensemble. 
In fact, the introduction of wave packet splitting into AMD was first motivated for 
the purpose to get a proper statistical properties with the fermionic caloric curve 
£;* = aT^ at low temperaturelll§'lll2l From a different point of view, Ohnishi 
and Randrup also introduced stochastic terms into molecular dynamics for quantum 
statisticsllHf The caloric curves in the region of liquid gas phase transition was 
calculated by SugawaandHoriuchi by employing AMD with an implementation of 
wave packet splitting]!^ 'llMt Fermionic molecular dynamics was also applied to the 
caloric curves!!^ The caloric curves obtained in these studies were drawn under 
the condition of fixed volumes or for a system under a confining potential. 

More recently, in Refs. 155), 156), Furuta and Ono performed AMD calculation 
to obtain constant-pressure caloric curves in which phase transition should be clearly 
identified. Wave packet splitting was considered with a density-dependent coherence 
time t(p) in this study. For a calculated microcanonical ensemble of a given energy 
E and a volume V, the temperature T is defined by using the kinetic energies of 
gas-like nucleons. The pressure P is obtained from the information of the reflections 
of particles at the boundary of the container. Then caloric curves are drawn as T{E) 
for different values of P. The result for the system of = 18 and Z = 18 is shown 



in Fig. [25) Nuclear liquid-gas phase transition, as a first order phase transition, is 
clearly seen in this result as the back-bending of caloric curves. 

The construction of microcanonical ensembles can also be regarded as a way to 
explore the many-body states at various excitation energies, as illustrated in Fig. 
S6| Above the ground state and low-lying excited states of the nucleus, there should 



be quite a lot of states which are characterized by the density of states W{E). The 
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Fig. 25. Constant pressure caloric curves calculated with AMD for the A = 36 system. The cross 
symbols indicate the equilibrium state corresponding to the reaction system from t = 80 to 300 
fm/c for ^"^CaH- ^"^Ca central collisions. The figure is taken from Ref. 156). 
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Fig. 26. Excited states of a many-nucleon system with A/" = Z = 18 confined in a virtual spherical 
container with a radius of 9 fm. The density distribution at each excitation energy (E* /A — 4, 10 
and 28 MeV) shows a snapshot taken from the AMD equilibrium calculation of Ref. 155). 



density profile at each of the excitation energies — 4, 10 and 28 MeV is a sample 

taken from the calculated microcanonical ensemble 11^ At low excitation energies, a 
single large nucleus is usually observed. The energy E"" /A — 10 MeV is in the region 
of liquid-gas phase transition where W{E) shows an anomalous behavior and each 
density profile typically shows several nuclei being related to multifragmentation. It 
should be noted that the system confined in a container does not become a gas of 
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nucleons even though the excitation energy is higher than the binding energy of the 
nucleus. When the energy is further raised to 28 MeV/nucleon, for example, the 
state may be regarded as a gas state but the gas is composed of clusters as well as 
nucleons. 

Thus it is now a great advantage of AMD that it can describe both dynamical 
reactions in heavy-ion collisions and virtually equilibrated systems reasonably well. 



In Ref. 156), we investigated the question whether equilibrium is really relevant in 
multifragmentation by comparing the details of reaction calculations and equilibrium 
calculations performed by the same AMD model. The calculations show that there 
exists an equilibrium ensemble which well reproduces the reaction ensemble at each 
reaction time t for the investigated period 80 < t < 300 fm/c in ^^Ca + ^^Ca cen- 
tral collisions at 35 MeV/nucleon, as far as fragment observables (fragment yields 
and excitation energies) are concerned. Thus the corresponding temperature and 
excitation energy (or the volume and pressure) can be identified at each reaction 



time. In Fig. 25, the path of the reaction form t = 80 to 300 fm/c is drawn by the 
cross symbols on the caloric curve figure. It is also important to note that there are 
some other observables which show discrepancies between the reaction and equilib- 
rium ensembles. ^^^^ These may be interpreted as dynamical effects in the reaction. 
In particular, the usual static equilibrium at each instant is not realized since any 
equilibrium ensemble with the same volume as that of the reaction system cannot 
reproduce the fragment observables. 

4.6. Symmetry energy effects in heavy-ion collisions 

The AMD simulations for heavy-ion collisions are useful not only to explain 
the experimental data but also to know what kind of information is reflected in the 
fragment formation. In particular, as demonstrated in the previous subsection, we 
may expect that statistical properties such as the equation of state can be extracted 
from the fragment observables even in dynamical collisions. 

In heavy-ion collisions with unbalanced neutron and proton numbers, the dif- 
ference between the neutron and proton motions is an interesting new degrees of 
freedom. The difference of flow pattern between neutrons and protons has been 
predicted to be sensitive to the density dependence of symmetry energy. From the 
viewpoint of liquid-gas phase transition, neutron-rich systems are quite interesting 
because new characters as two-component systems are expected. Namely, the gas 
part of the system is more neutron-rich than the liquid part, which can be called 
fractionation or distillation. This effect of isospin fractionation/dist illation should 
be observable in the neutron-to-proton ratio of produced fragments. 



In Refs. 158), 159), the fragment yields were analyzed in the AMD simula- 
tions for multifragmentation reactions of the central collisions of Ca isotopes at 
35 MeV/nucleon, in order to see how the fragment isospin composition is related 
to the symmetry energy term of the effective interaction adopted in the calculation. 
For the fragment yields Yi{N, Z) in the reaction z at t = 300 fm/c, it was found that 
the isoscaling relation 

Yj{N, Z)/Yi{N, Z) oc e^^^-^+^^^-^ (4-6) 
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is satisfied for any two reaction systems i and j which are different in the proton- 
to-neutron ratios. Isoscahng is expecte d u nder an equihbrium assumption and has 
been observed in the experimental datall^Zl Isoscahng is equivalent to the statement 
that the fragment yields are expressed as 

y,(7V, Z) = exp[-K(N, Z) + a,N + ^.Z + 7,] (4-7) 

by using a function K{N, Z) that is independent of the reaction system i. If the 
equilibrium is relevant to the reaction, K(N^ Z) should contain a term (Csym/r)(A^ — 
Z)V(7V + Z), where Csym is a kind of symmetry energy which may depend on A or 
Z in principle. Based on this assumption, the isoscahng parameter is related to the 
symmetry energy by 



4Csym 



{Z/MZ)f-{Z/A,{Z)f , (4-8) 



where Z/Ai{Z) represents the mean isospin asymmetry of the fragments for each 
given Z in the reaction i. 




Fig. 27. Density dependence of the symmetry energy of nuclear matter for the Gogny force (sohd 
hne) and for the Gogny- AS force (dashed hne). 



AMD calculations were carried out with the Gogny force and the Gogny-AS force 
which are different in the density dependence of the symmetry energy as shown in 
Fig. [27| The result shows that the isoscahng parameter a actually depends on the 

158) Furthermore, 



28 



density dependence of the symmetry energy as shown in Fig. 
the value of Cgym/T extracted from the simulation result is found to be almost 
independent of the fragment size A or z'J^ which suggests that bulk properties 
are reflected in the fragment isotope yields rather than the symmetry energy for the 
ground state binding energies that depends on A due to the surface effect. To explain 
the obtained result, Cgym is identified with the bulk symmetry energy at about |po, 
with po being the saturation density of nuclear matter, and the temperature should 
be T ^ 3.4 MeV. These values of the density and the temperature are reasonable as 
the condition for fragmentation. 
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Fig. 28. The fragment yield ratio between the AMD simulations of central ^°Ca + ^°Ca and 
^°Ca + ^°Ca collisions at 35 MeV/nucleon, at time t — 300 fm/c. The top and bottom panels 
show, respectively, the results obtained using the Gogny and Gogny-AS forces. The extracted 
isoscaling parameters are a — 1.82 it 0.06 and (3 — —2.23 ± 0.08 for the Gogny force, and 
a = 1.64 ± 0.05 and /3 = -2.09 =b 0.07 for the Gogny-AS force. The figure is taken from Ref. 
158). 



It should be noted that the above analysis has been done for the fragment yields 
at t = 300 fm/c. In order to compare with experimental data, the effects of the 
decay of these primary fragments should be carefully considered. 

§5. Summary and perspective 

In nuclear systems, cluster aspect is one of the essential features as well as 
mean- field aspect. Coexistence of cluster and mean- field aspects brings a variety of 
phenomena as functions of excitation energy and isospin degrees of freedom. For 
usual theoretical models, it is not easy to describe both behaviors of independent 
single nucleons in a mean field and spatially correlating nucleons in clusters. The 
AMD method is one of the theoretical approaches that can describe those two kinds 
of nature. In AMD, single-particle wave functions are written by localized Gaussian 
wave packets whose dynamics expresses assembling and disassembling of nucleons. 
The method has been applied to investigate nuclear reactions and structures and it 
has been proved to successfully describe a variety of phenomena in general nuclei. 

In this paper, we reviewed the AMD approach and its application to nuclear 
structure and reaction studies. To show applicability of AMD some topics studied 
with time-independent and time-dependent versions of AMD were explained. In 
the applications of time-independent AMD to nuclear structure studies, structures 



Antisymmetrized molecular dynamics 



43 



of neutron-rich nuclei such as Be, C, F, Ne, and Mg isotopes were described fo- 
cusing on cluster aspects. The results suggested a variety of structures appear in 
unstable nuclei as well as in stable nuclei. Deformation and cluster phenomena in 
Z N nuclei in p- and sd- shell regions were also described. The applications 
of time-dependent AMD contain various topics such as fragmentation in heavy-ion 
collisions, nuclear responses and virtual systems in thermal equilibrium. AMD cal- 
culations successfully described multifragmentation which is one of the remarkable 
phenomena in heavy-ion collisions. Comparisons with predictions by other models 
and/or experimental data suggest the important balance between the single-particle 
motions and the many-body correlations to form clusters and fragments in these 
phenomena. The AMD approach is suited to link the reaction observables to the 
equilibrium properties of nuclear matter such as liquid-gas phase transition and the 
equation of state of asymmetric nuclear matter. 

Success of those studies using the AMD approach greatly owes to characteristics 
of the AMD model, for instance, advantages hsted below. 

• It is able to describe cluster and mean-field aspects without assuming existence 
of clusters nor mean fields. 

• It is applicable to both static and dynamical problems. 

• It it applicable to general nuclei with given proton and neutron numbers. 

• Center-of-mass motion can be exactly extracted. 

In addition to the above advantages, the model can be easily extended because of 
flexibility of AMD wave functions. For instance, parity and angular-momentum pro- 
jections and superposition of wave functions are performed in structure calculations, 
and stochastic collisions are incorporated in calculations of heavy-ion collisions. 

Present structure studies with the AMD method cover light-mass regions of 
nuclear chart up to p/-shell nuclei. It is a future problem to apply the method to 
further heavy-mass regions and progress systematic studies covering wide regions of 
the nuclear chart. 

Origins of cluster formation and breaking should be clarified from the point of 
view of nuclear force. Unfortunately, the present AMD framework is model calcu- 
lation and it requires phenomenological effective nuclear interactions. To calculate 
nuclear systems based on realistic nuclear forces is one of the important issues in 
nuclear physics. In fact, there are many attempts of ah initio calculations though 
practical ah initio calculations are still limited in very light systems. For system- 
atic study covering a wide mass number region, model calculations are efficient. To 
achieve model calculations starting from realistic nuclear forces a main problem is 
how to deal with complicated many-body correlatio ns. One of th e promising methods 
is the unitary correlation operator method (UC0M)ll£2}'ll£ll'[l£l} recently proposed to 
incorporate efficiently short-range and tensor correlations in structure models such 
as FMD. To understand nuclear structure from fundamental point of view, more 
sophisticated versions of the AMD method would be required. 

Another issue of the AMD approach would be wide application to reaction phe- 
nomena. Although the AMD method has been extensively applied to violent reac- 
tions above ten MeV/nucleon, its applications to lower-energy reaction phenomena 
are limited. Motivated by recent progress of experimental studies, low-energy re- 
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actions such as resonances, fusion/capture, and transfer reactions are interesting 
problems to be solved concerning cluster phenomena in unstable nuclei. 
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